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ABSTRACT 

o 

' Context. Ray-like features observed by coronagraphs in the wake of Coronal Mass Ejections (CMEs) are sometimes 

interpreted as the white light counterparts of current sheets (CSs) produced by the eruption. The 3D geometry of these 
■ ray-like features is largely unknown and its knowledge should clarify their association to the CS and place constraints 

on CME physics and coronal conditions. 

Aims. If these rays are related to field relaxation behind CMEs, therefore representing current sheets, then they should 
be aligned to the CME axis. With this study we test these important implications for the first time. 
Methods. An example of such a post-CME ray was observed by various coronagraphs, including these of the Sun Earth 
Connection Coronal and Heliospheric investigation (SECCHI) onboard the Solar Terrestrial Relations Observatory 
(STEREO) twin spacecraft and the Large Angle Spectrometric Coronagraph (LASCO) onboard the Solar and 
Heliospheric Observatory (SOHO). The ray was observed in the aftermath of a CME which occurred on 9 April 
2008. The twin STEREO spacecraft were separated by about 48° on that day. This significant separation combined 
i-C ■ with a third "eye" view supplied by LASCO allow for a truly multi-viewpoint observation of the ray and of the CME. 

We applied 3D forward geometrical modeling to the CME and to the ray as simultaneously viewed by SECCHI-A and 
B and by SECCHI-A and LASCO, respectively. 

Results. We found that the ray can be approximated by a rectangular slab, nearly aligned with the CME axis, and 
much smaller than the CME in both terms of thickness and depth (~ 0.05 and 0.15 Rq respectively). The ray electron 
density and temperature were substantially higher than their values in the ambient corona. We found that the ray and 
CME are significantly displaced from the associated post-CME flaring loops. 

Conclusions. The properties and location of the ray are fully consistent with the expectations of the standard CME 
theories for post-CME current sheets. Therefore, our multi-viewpoint observations supply strong evidence that the 
^ ' observed post-CME ray is indeed related to a post-CME current sheet. 
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^ , 1. Introduction existence is in the form of White Light (WL) rays seen by 

__l ■ coronagraphs behind CMEs (e.g., Ko et al. 2003; Webb et 

■ CMEs result from a major restructuring of the corona. al 2003; Vrgnak et al 2009) The rays are formed in the 

^ ; During such events, magnetic fields and plasma are ex- wake of CMEgj they extend tQ geveral solar radii and per _ 
j. . pelled from the Sun causing dramatic changes m the coronal sigt from seyeral hours up to few dayg The detection of 
m ;> ■ medium both within and around the erupting volume. yery hot (> 3 MK) plasmas in several of these WL rays by 
k> Essentially all CME models predict the formation of the TJltraViolet Coronagraph Spectrometer (UVCS; Kohl 
r ^ ■ a post-CME current sheet (CS) underneath the erupting et aL 1995) further sup p 0rte d the possibility that such 
C$ ; magnetic flux rope (e.g. see also the review of Forbes and structures represent CSs (e.g., Ko et al. 2003; Raymond 
Lin 2000). As the magnetic flux rope U ascends, it pushes et al 2 003; Bemporad et al. 2006; Lee et al. 2006; Lin et 
aside overlying magnetic fields. These fields will pinch off al 2007 ) Coronal CSs are places where one could expect 
underneath the erupting rope forming a neutral CS, i.e. a higher heating rates, and therefore higher temperatures, 
surface separating fields of opposite magnetic polarity. The Moreover, such rays are not quiescent but they often ex- 
faster the eruption proceeds, the higher the rate at which hibit dynamics such as propagating blobs through the ray 
new field lines accumulate at the CS. Such post-eruption ( eg ^ S heeley and Wang 2007). Such dynamic phenomena 
CSs should be distinguished from the CSs which could exist could result from re connections taking place in a CS. A re- 
fee/ore the eruption in some CME models. cent review of cs observations made by LASCO and UVCS 
There exist several pieces of observational evidence sup- i s g i ven m Vrsnak et al. (2009). An analytical Petschek-like 
porting the existence of post-CME CSs both in the inner curre nt sheet model was found in a good agreement with 
and in the outer corona. A major "smoking gun" for their tne observations. Closer to the solar surface, observations of 
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cusp-like hot structures and descending voids above erup- 

1 Flux ropes are currently widely recognized as the major con- tion sites and arca des provided further support for the exis- 
stituent of CMEs. 



2 



Patsourakos and Vourlidas: STEREO SOHO observations of a post-CME current sheet proxy 



tence of post-CME CSs (e.g., McKenzie and Hudson 1999; 
Innes et al. 2003; Asai et al. 2004; Sheeley et al. 2004). 

However there are a couple of reasons which argue 
against the possibility the observed rays could represent 
CSs. First, standard plasma theory predicts that CSs 
should be very narrow (a few meters or less) for typical 
coronal conditions (e.g., Litvinenko 1996; Wood & Neukirch 
2005). However, the coronagraph observations suggest CS 
widths at a fraction of the solar radius 0.2 — 0.8i? Q : Vrsnak 
et al. 2009). One way to overcome this apparent contradic- 
tion is to invoke anomalous processes (e.g., turbulence, hy- 
peresistivity) which could substantially broaden CSs (e.g., 
Lin et al. 2007; Bemporad 2008). Another possibility is that 
the macroscopic CS could be in reality the superposition of 
several small hot-spots, patches where small-scale reconnec- 
tions take place (e.g., Klimchuk 1996; Cargill et al. 2006; 
Linton and Longcope 2006; Onofri, Isliker, & Vlahos, L. 
2006, Vrsnak et al. 2009). 

A second difficulty comes from the very nature of all 
past ray observations which employed a single view point. 
This does not allow to determine the full 3D geometry (i.e., 
width depth and height) of the rays and their exact rela- 
tionship with the associated CMEs. All we could measure 
are quantities projected in the plane of sky (POS) of the 
corresponding instrument. Even more importantly one has 
to worry that projection effects in the optically thin corona 
could cause one to falsely believe that the observed rays rep- 
resent post-eruption CSs. This is particularly true during 
CMEs which cause reconfigurations of significant coronal 
volumes in terms of rotations, deflections, displacements 
etc of structures both within and around the CME. For in- 
stance the outer CME legs or structures around the erupt- 
ing volume, could be falsely identified as candidate CSs 
when observed at a favorable angle; they may appear to 
trail the CME in projection as viewed from a single view- 
point. However, such structures may be far from the erup- 
tion core where the CS should lie. 

Therefore we do not know the "true" width and depth 
of these rays, not to mention their shapes. Moreover, we 
currently ignore how they are arranged in 3D space with 
respect to the associated CME and whether they represent 
a projection effect. From the standard theory for the forma- 
tion of post-CME CSs one would expect them to be quasi- 
aligned with the CME and to be bounded within it (i.e., 
to occupy a smaller volume than the CME itself). Clearly 
observations of post-CME rays from multiple and distinct 
viewpoints are required in order to address the above im- 
portant questions. 

The launch of the STEREO mission in late 2006 ini- 
tiated an avenue for the study of the solar corona in 3D. 
The mission consists of two almost identical spacecraft, one 
orbiting ahead (STA) and the other orbiting behind (STB) 
Earth's orbit. The two spacecraft drift apart at a rate of 
ps 45° per year. With this paper we supply the first de- 
tailed 3D analysis of a post-CME ray using STEREO and 
LASCO data. Our analysis supplies some new and impor- 
tant geometrical and physical evidence that the observed 
ray could represent a CS. 

Our paper is structured as follows. An overview of the 
CME-ray observations is given in Section [5] Section [3] de- 
scribes the 3D geometrical modeling applied to the CME 
and to the ray whereas in Section U we determine various 
physical parameters (density and temperature) of the ray. 



Finally, Section [5] summarizes our findings and discusses 
their implications. 

2. Overview of the Observations 

We analyze here a CME event which took place on 9 April 
2008. The event was initiated in the active region (AR) 1098 
of the NOAA classification. There was no GOES flare asso- 
ciated with this CME: this could be due to the fact that the 
source active region was partially occulted behind the West 
limb as seen from the Earth. The event was observed by 
the Solar and Heliospheric Observatory (SOHO), Hinode, 
the Transition Region and Corona Explorer (TRACE) and 
STEREO. A synopsis of these observations was presented 
by Reeves at al. (2008), Landi et al. (2010) and Savage et 
al. (2010). 

We have posted a movie (moviel.mpg) showing the low 
coronal development of the event. This movie consists of 
composite EUVI-195 A channel and COR1 Extreme Ultra- 
violet and WL images respectively taken by STA; a cou- 
ple of representative snapshots from this movie is given in 
Figured] the FESTIVAL software of Auchere et al. (2008) 
was used to generate the composite images. 

The event started around 09:10 UT with the activation 
and the lift off of a prominence. The eruption moved ini- 
tially parallel to the surface towards the southeast and was 
diverted outwards along a southwestern path. Such deflec- 
tion may have been caused by the large equatorial coronal 
hole south of the host active region. The strongly non-radial 
propagation of the early CME can also be evidenced in the 
upper panel of Figure[T]where we note a significant displace- 
ment between the axis of the emerging CME into the COR1 
field of view (FOV) and the position angle of the host AR. 
The erupting prominence eventually gave rise to a typical 
3-part CME (e.g., upper panel of Figure [lj . A rather faint 
linear feature emitting in soft Xrays in the trail of the erupt- 
ing CME was observed by the X-ray Telecope (XRT) on 
Hinode (Reeves et al. 2008; Savage et al. 2010). Combined 
analysis of the XRT and the EUV Imaging Telescope (EIS) 
Hinode observations showed that this feature had temper- 
atures in the range « 6-10 MK (Landi et al. 2010). This 
ultra-hot linear structure in the wake of a CME may have 
been associated with a post-CME CS. Indeed, XRT count- 
rates as predicted from a Petchek-like current sheet model 
were found in a good agreement with the weak signal of 
the mentioned above current-sheet like feature (Ko et al. 
2010). Once the CME exited the COR1 FOV a WL ray 
was formed at its wake (e.g., lower panel of Figure [TJ. Also 
Savage et al. (2010) observed with XRT mass flows (e.g., 
outflows, supra-arcade downflows) around the mentioned 
above faint linear feature which supplied evidence that this 
structure possibly represented a post-CME CS. 

We concentrate here on WL observations of the corona 
taken by the COR2 coronagraphs (FOV from 2.5-15 i? Q ) 
of the Sun-Earth Connection Coronal and Heliospheric 
Investigation (SECCHI; Howard et al. 2008) instrument 
suite on-board the twin STEREO spacecraft and by the C2 
coronagraph (FOV from 2.5-6 Rq) of the LASCO experi- 
ment on SOHO (Brueckner et al. 1995). On the day of our 
observations the twin STEREO spacecraft were separated 
by w 48°. Our combined multi- viewpoint observations al- 
lowed to monitor and characterize in 3D the evolution and 
the structure of the CME and of the post-CME ray. Raw 
data were first treated with the standard data reduction 
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routines. COR2 obtains both total WL brightness (TB) 
and polarized WL brightness (PB) images, whereas C2 nor- 
mally obtains only TB images. Typical image cadence is of 
the order of 10-20 minutes. There were no UVCS observa- 
tions of the ray as the instrument was observing lower in 
the corona (Landi et al. 2010). 

Figure [5] contains representative TB images through- 
out the event for COR2 on STA (COR2A), C2 and COR2 
on STB (COR2B) respectively; the full sequences can be 
found in the corresponding online movies (movie2.mpg, 
movie3.mpg and movie4.mpg). From each image in these 
figures and the associated movies a pre-event TB image was 
first subtracted creating a base-difference image sequence. 
Therefore, bright (dark) areas correspond to mass increase 
(decrease) respectively. 

Several remarks can now be made. The CME has al- 
ready emerged into the COR2A and C2 FOV around 11:00; 
this happens later for COR2B and hints to a CME being 
closer to the STA and SOHO plane-of-sky (POS). Once it 
appears in the coronagraph FOVs the CME has the canon- 
ical 3-part CME structure with a bright front, dark cavity 
and bright core, i.e. it is a typical flux-rope CME. The 
CME core exhibits significant fine structure and consists of 
many threads as can be seen better in the C2 images. The 
CME seems to undergo a rotation as it expands. Indeed, 
Thompson, Torok & Kliem (2010) studied the rotation 
of the erupting prominence that was associated with our 
CME. Using triangulations of several threads of the erupt- 
ing prominence observed by STA and STB they determined 
its rotation profile for heliocentric distances « 1.5-3.3 R@. 
They found that the prominence rotated by « 120° rela- 
tive to its preeruption position. The CME starts to exit 
the SOHO FOV around 13:00 UT and the STA and STB 
FOVs around 15:52 UT and 16:52 UT, respectively. While 
this happens we see a WL ray forming in the wake of the 
CME in C2 (e.g., from 14:50 UT onwards in Figure© and 
in COR2A (e.g., from 15:52 UT onwards in Figure©; there 
is no evidence for such a ray in STB though. The WL ray 
seems to undergo a rotation as the CME itself is rotating 
during its expansion. The WL ray connects to the bottom 
part of the CME cavity and seems to be formed when lin- 
ear features joining the CME core start to pinch off as the 
CME propagates outwards. The ray was visible for at least 
one day after its formation. 

From the above we conclude that the observed WL ray 
could be a candidate CS. One needs to always keep in mind 
that strictly speaking a CS is a surface separating magnetic 
fields of opposite direction: this is something which can not 
be determined by our observations or by any other similar 
observation given the lack of direct magnetic field obser- 
vations in the corona. The 3D analysis of this paper will 
supply rather strong indirect evidence that the observed 
WL ray could have been a post-CME CS. 

A first set of clues about the line of sight (LOS) extent of 
the post-CME ray and of the CME could be drawn from the 
analysis of the COR2 TB and PB images shown in Figure 
[3J The theory of Thompson scattering (Billings 1966) which 
describes the formation of the WL emissions observed by 
coronagraphs predicts that PB decreases very fast when the 
emitting electrons lie away (i.e., w ±20°) from the POS of a 
given instrument (Vourlidas and Howard 2006). Therefore 
PB tends to emphasize structures close to a given POS; 
this is not the case for TB which has a more shallow fall-off 
with distance from the POS. 



What we note from Figure [3] is that the ray is absent 
in the STB PB and TB images whereas it can be seen in 
STA PB and TB. This means that the ray should lie close 
to the POS of STA. 

Moreover, the LOS extent (depth) of the ray cannot be 
too large (i.e., comparable to the CME width). Otherwise, 
the ray would be visible in STB in cither TB or PB since 
the CME itself was seen by both spacecraft in TB (Figure 
©) and in PB (Figure G2 although less pronounced in STB). 
The above observations supply a rather strong indication 
that the ray could not have been due to a projection effect 
(e.g, widely separated CME legs lining up along the LOS) 
but it is rather a structure with limited extent along the 
LOS. This is further demonstrated and quantified in the 
next Section. 

3. 3D Geometrical Modeling of the CME and of 
the post-CME ray 

In order to determine the 3D structure of the post-CME 
ray and to infer its 3D position with respect to the CME 
we performed 3D geometrical modeling of both the CME 
and the ray. 

We first start with the CME modeling. For this task we 
used the 3D forward model of Thernisien et al. (2006, 2009) 
(rtsccguicloud.pro and rtraytracewes .pro ssw rou- 
tines) . The model has been successful in fitting the 3D mor- 
phology of CMEs (Thernisien et al. 2006, 2009) observed 
either from one viewpoint (LASCO) or from two viewpoints 
(STA and STB). The latter approach deemed particularly 
powerful in strongly constraining the CME parameters. We 
used TB images from COR2 on STA and STB. 

Here, we model the observed CME with a spherical shell 
attached to two conical legs, a so-called "croissant" model. 
The spherical shell is meant to reproduce the upper part of 
the CME, whereas the conical legs are an approximation for 
the CME's lower sections, underneath the spherical shell. 
Such a model emulates the outline of a flux-rope, widely 
accepted as the magnetic structure of CMEs. The model is 
formulated in terms of the following free parameters: po- 
sition (longitude (9), latitude (</>) of the legs on the solar 
surface ; tilt (7) of the model baseline with respect to the 
equator; height above the solar surface and aspect ratio 
(=height/radius) of the spherical shell. The model is ra- 
dial, i.e. it lies vertical to the solar surface. A schematic of 
the model along with its free parameters can be found in 
Figure |4] 

The model parameters were varied until we achieved a 
visually satisfactory agreement between the observed CME 
(in TB) and the wireframe projection of the model was 
achieved simultaneously for STA and STB (see left and 
right columns of Figure [5] respectively). We note that the 
model reproduces fairly well both the body and the legs of 
the CME. We plot in the middle column of Figure [5] syn- 
thetic TB images from the CME model for STA and STB 
to illustrate the model fitting further . A uniform density 
of 10 6 cm~ 3 was assumed; its precise value does not affect 
the purely geometrical analysis of this Section. Nevertheless 
the deduced geometrical information here will be used to 
determine the ray density in the following Section. The fig- 
ure shows that the model appears to reproduce some of the 
internal structure of the CME. The WL renderings from 
the model show some fine structure within the CME which 
is also present in the actual observations. 
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To estimate uncertainties in the model parameters we 
perturbed each parameter of the "best-fit" model shown in 
Figure [5] by keeping the remaining parameters frozen. Each 
parameter was varied until a visually unacceptable solu- 
tion (i.e., projections) was found; this supplied a measure 
of the uncertainty for each parameter. The fitting proce- 
dure was purely visual. Namely, we derived a set of pa- 
rameters that give the best visual correspondance between 
the model and the two- viewpoint observations. The "best- 
fit" parameters along with their uncertainties of the CME 
model are given in Table 1. We note that the uncertainties 
are quite small, which underlines the importance of using 
multi-view point data to tightly constrain the CME model. 
This is particularly true for the CME longitude which indi- 
cates the direction along which the CME propagates. Given 
the Carrington longitude of STA was ~ 121° on this day, 
the deduced Carrington longitude of the CME (w 194°) 
implies that the CME was about 17° away from the STA 
POS and therefore almost 11° and 67° degrees from the 
POS of SOHO and STB respectively. Therefore the CME 
was mainly an Earth and STA event. A schematic with the 
orientation and width of the fitted CME along with the 
locations of STB, Earth and STA and their corresponding 
POS can be found in Figure [SI 

We then performed 3D geometrical modeling of the 
COR2A and C2 ray observations. We used again the 
Thernisien et al. (2006, 2009) approach but with a different 
3D geometrical model. A rectangular slab was used instead. 
This is a reasonable approximation for the ray shape since 
in both viewpoints it looks like an elongated with an almost 
constant cross-sectional area surface. The standard theory 
of solar eruptions (e.g., Forbes 2000) predicts post-CME 
current sheet shapes resembling rectangular slabs. Our ray 
model was formulated in terms of the longitude and lat- 
itude of the slab axis, and the slab width and depth. A 
schematic of the model along with its free parameters can 
be found in Figure |H Supposing for the moment that the 
ray represents a CS then its depth would correspond to 
the longitudinal extent of the post-CME arcade (along the 
erupting neutral line) whereas the ray width would corre- 
spond to the thickness of the current layer where opposing 
magnetic fields from both sides of the neutral line come 
into a diffusion region. For a ray viewed almost edge-on its 
depth will be comparable to the LOS extent. Since the ray 
width and depth are similar, the slab has a nearly square 
projection so the tilt becomes ambiguous (i.e., various tilts 
can fit the data equally well). We therefore "fixed" the ray 
tilt to the "best-fit" tilt of the CME and derived a width 
and depth shown in Table 1. 

Figure [7] contains our "best-fit" model of the ray. As 
done for the CME the "best-fit" model was determined 
by visual comparison between the model predictions and 
the actual multi-viewpoint observations. As can be seen 
by this Figure our model reproduces the ray orientation 
and projected size in both viewpoints. The uncertainties 
were derived the same way as done for the CME. The ray 
parameters along their uncertainties are given in Table 1. 
Again, the simultaneous fitting of the two viewpoints sup- 
plied strong constraints on the model parameters with the 
exception of the ray tilt angle due to its symmetry. For in- 
stance we note that increasing the ray depth and width by 
a factor of two we get unacceptable solutions (Figure [8J ; 
the ray is too wide in the SOHO view. Two important re- 
marks can be now made. First, the ray has a rather small 



Table 1. CME and ray fitting results. 



CME CS 



Carrington longitude(°) 
Carrington latitude(°) 


194 ± 5 


197 ± 5 


-17 ± 3 


-20 ± 3 


tilt(°) 


4 ± 7 


4 ± 7 


aspect ratio 


0.16 ± 0.03 


N/A 


angular width(°) 


18 ± 4 


N/A 


depth (R Q ) 


N/A 


0.15 ±0.08 


thickness (Rq) 


N/A 


0.05 ±0.03 


density (cm~ 3 ) 


N/A 


3xl0 7 ±5.3 x 10 6 



width (0.05 Rq) and depth (0.15 Rq) which represent a 
small fraction of the solar radius. This can be contrasted 
with the large radius of the CME (e.g., around 1.8 Rq at 
14:52 UT). Second, the location of the ray (i.e, longitude 
and latitude) is not very far off from the location of the 
CME. For instance the CS central longitude differs by only 
3° from that of the CME. This small difference is almost 
within the error bars and certainly much smaller than the 
CME width of w 18° This means that the ray is almost 
perfectly aligned with the CME and lies close to its center. 
Indeed the ratio between the CME to ray volumes is « 50, 
meaning that the ray represented only a small fraction of 
the CME volume. We will discuss the important implica- 
tions of the above findings in Section [5] 

4. Determination of the Ray Density Structure 

Since the WL intensities in coronagraph images are pro- 
portional to the LOS integral of the electron density (e.g., 
Billings 1966) we used our observations to infer the den- 
sity structure of the ray. The common assumption made in 
such determinations (e.g., Poland et al. 1981; Vourlidas et 
al. 2000; Vrsnak et al., 2009) is that the structure of inter- 
est lies in the POS of the corresponding instrument. Then, 
local densities are inferred by assuming a LOS extent; a 
common assumption made is that it is proportional to the 
projected width of the structure. 

Our measurements of the 3D shape and location of the 
ray allowed to relax the above two assumptions by supply- 
ing the distance of the ray from the STA POS (« 15°). 
Since the ray orientation is degenerate the LOS extent of 
the ray can be in principle anywhere between the calcu- 
lated ray thickness (0.05 Rq) and depth (0.15 Rq). We 
make the assumption it is 0.10 i?©, i.e. the median value of 
the above interval. Therefore, the LOS extent uncertainty 
is ±O.O5i?0, which translates into 50% density uncertainty. 
In such a way we were able to determine "true" densities. 

The determination of the ray densities involved the fol- 
lowing steps. First, we created an excess mass image of the 
STA ray for the TB observation 18:07:54 UT UT by sub- 
tracting a pre-event image (taken at 09:22:54 UT). Both im- 
ages were fully calibrated. The conversion from Mean Solar 
Brightness to electrons/cm -2 was done following the stan- 
dard assumptions of Thompson scattering (e.g., Vourlidas 
et al. (2000) for details) and the distance of the ray from 
the STA POS. 

The excess mass image is shown in Figure [5] White 
(black) regions represent places with excess (deficit) mass 
with respect to the pre-event reference image. Several 
points along the ray were then manually selected (dashed 
line in Figure ■ To enhance the ray signal the mass mea- 
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surements along the ray were binned every 0.02 Rq. Then 
the excess mass of the ray was converted into a density by 
using its LOS extent from our 3D modeling. We note here 
that the derived densities are likely lower limits to the true 
densities. If the CME removes significant part of the over- 
lying streamer, then the pre-event subtraction will result in 
lower excess brightness than expected. This is insignificant 
in our case, since there is only a diffuse preexisting struc- 
ture at the location of our ray and the observed depletions 
are minimal at the ray location (especially on the north- 
ward side). Uncertainties due to photon counting statistics 
can be safely ignored given the very large count rates we 
have in WL at the heights of our observations (typically 
above thousand counts per pixel). The uncertainty in the 
excess mass measurements of the COR2 instruments is es- 
timated to be 20 %. The total uncertainty arises from the 
quadratic combination of these uncertainties (mass deter- 
mination and LOS depth) under the assumption they are 
independent (i.e., we took the average of the quadratic sum 
of the individual uncertainties) . 

The resulting density profile along the ray is given in 
Figure [TDJ Note also that the distances along the ray are 
also "true" in the sense we used our 3D CS modeling to 
convert POS distances to heliocentric distances. We plot 
densities only in the range 3.2-3.8 Rq because the ray be- 
comes too diffuse at larger distances to be discerned from its 
surroundings. Several remarks can now be made. First, the 
derived density values are more than two orders of magni- 
tude larger than the densities of the pre-CME corona at the 
heights of our observations. This can be seen in Figure [10] 
where we plot densities from the inversion of the LAS CO 
C2 PB sequence taken at 10:06 on the day of the event. 
The sequence was taken before the emergence of our CME 
into the C2 FOV and gives an estimate of the densities 
in the pre-CME corona. The pre-CME densities were cal- 
culated at the same position angle as the ray densities in 
Figured 1101 The deduced ray densities are within the range 
» 10 7 -10 8 cm determined for several WL rays by Vrsnak 
et al. (2009). Applying the same procedure to COR1 data, 
we found a ray density of pa 3 x 10 8 cm -3 at a distance of 
about 1.6 Rq. Again the ray emission becomes too faint at 
larger heights to be reliably measured. Second, we note that 
the ray density does not substantially decrease with height 
along its length (only a factor w 2 over a distance of ~ 0.5 
Rq). This implies reduced scales heights and therefore ele- 
vated temperatures within the ray. Indeed, we can see from 
Figure [TU] that the ray density drop-off can be fairly well 
reproduced by a hydrostatic density profile at a uniform 
temperature of 2 MK. In doing this we used a profile of the 
following form : 



n{h) — rioexp(—h/H), (1) 

where h is the ray height above the solar surface, no the den- 
sity at the base of the ray and H the isothermal scale-height 
(i.e., oc to temperature). Empirical models of the quiet 
corona indicate temperatures of < 1 MK at the heights of 
our observations (e.g., Withbroe 1988). We note here that 
the temperature is inferred under the assumption of a hy- 
drostatic equilibrium. Such an assumption may not be fully 
justified for post-CME rays which are often characterized 
by dynamic phenomena such as plasma blobs. 



5. Discussion and Conclusions 

This paper presents the first multi-viewpoint observations 
of a post-CME ray. Utilizing multiple view points together 
with detailed 3D geometrical modeling supplied new im- 
portant insight into the 3D structure of a post-CME ray 
and its relationship with the associated CME. Our main 
findings are the following: 

— The ray can be fitted by a rectangular slab with small 
(w 0.15i? Q ) depth and (« 0-05.Ro) thickness 

— The ray axis is almost aligned with that of the CME 

— The ray lies within the CME volume and it occupies a 
much smaller volume than the CME 

— The ray is characterized by elevated (> 100) densities 
with respect to the ambient corona 

— The ray density drop-off with height is consistent with 
a high (w 2 MK) temperature, significantly hotter than 
the ambient corona at the same heights. 

Our ray density and temperature findings are fully in 
line with previous post-CME ray observations with UVCS 
and LASCO which showed they are characterized by ele- 
vated densities and temperatures. It is reasonable to expect 
such behavior within CSs. The pinch-off of field lines result- 
ing into a CS would compress the plasma thereby increasing 
its density. Also significant heating may occur in CSs since 
they are characterized by enhanced current densities. 

Undoubtedly the most important implications of our 
study arise from the 3D analysis of the event. We have de- 
termined the 3D shape of a post-CME WL ray and deduced 
its width and depth. Figure QTj contains an overlay of the 
CME and the WL ray central longitudes and latitudes as 
determined from the 3D modeling of Section 3 on a MDI 
magnetogram taken on 1 April 2008, i.e. close to the cen- 
tral meridian passage of the host AR as viewed from SOHO. 
We first note that both the CME (red cross) and WL ray 
(blue cross) are significantly misplaced towards the south- 
west direction with respect to the host AR. This suggests 
a non-radial evolution of the eruption in the inner corona 
and verifies the visual impressions of Section 2 (e.g., Figure 
[1]). Moreover, the 3D ray depth (equal to the vertical blue 
dashed lines of Figure ITT]) is comparable to the length of 
the neutral line of the host AR. If the ray corresponds to 
a CS, we would expect its depth to be comparable to the 
length of the erupting neutral line as seems to be the case. 
However, the deduced ray width (« 0.05i? Q ) is still many 
orders of magnitude bigger than the expectations of stan- 
dard plasma theory. 

We also found that the CME and the ray were al- 
most perfectly aligned. This supplies strong evidence for 
a tight and causal relationship between CME expansion 
and ray(CS) formation. The ray(CS) should trail the CME 
if it results from the pinch off of magnetic field behind 
the eruption. The fact that the WL ray could represent 
a post-CME CS could help removing the degeneracy in 
the ray orientation discussed in Section 3 since they ray 
should lie in a plane perpedicular to the plane of CME 
expansion. The fact that the ray and CME are not per- 
fectly aligned can be rather easily anticipated in terms of 
the various re-arrangements (rotations, deflections etc) tak- 
ing place during the event. Furthermore, according to the 
standard model for solar eruptions one expects that the 
post-eruption CS in the inner corona should lie above the 
post-flare loops which in turn are straddling the neutral 
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line of the erupting AR. The fact that the candidate CS 
(i.e., the WL ray) is misplaced with respect to the erupting 
neutral line suggests the possibility of a "broken" current 
sheet. This introduces a new, and possibly important, con- 
sideration for CME models. 

As a further visualization of the intimate relationship 
between the ray and the CME Figure [T3] presents WL ren- 
derings of the ray and the CME from our modeling taken 
from two different views: one from the ecliptic plane and 
another 34° above the ecliptic. The latter is a configura- 
tion which will be possible with the planned Solar Orbitcr 
mission. What we note in this Figure is that the ray is 
much smaller than the CME and it is located between its 
inner legs. It seems therefore likely according to the stan- 
dard model for solar eruptions that once the CME leaves 
the FOV its legs will pinch off and give rise to a ray (CS). 

In summarizing, our study clearly shows that the ray 
lies where it should be (aligned with the CME) and it is 
as big as it should be (depth similar to the erupting active 
region neutral line) if it had been a CS. Therefore, all our 
3D findings give strong support to the possibility that the 
observed post-CME ray is related to a post-CME CS. 
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Fig. 1. Composite of EUVI-A 195 (inner FOV) and COR1-A (outer FOV) showing the emerging CME (upper panel) 
and the resulting WL ray (lower panel). The temporal evolution is shown by moviel.mpg. 



8 



Patsourakos and Vourlidas: STEREO SOHO observations of a post-CME current sheet proxy 





\ ^ 


w 




V 




v \Jf| 








H 






^ VI 










KB 




m 




- 









Fig. 2. Representative total brightness (TB) snapshots of the CME and of the WL ray as observed by COR2A, LASCO 
C2 and C0R2B (upper, middle and lower panel respectively). Pre-event TB images taken at 01:27:13 UT (C2) and 
09:22:00 UT (COR2A-B) were subtracted from each C2, COR2A, and COR2B frame. The C2 field of view spans 2.2-6 
Rq and the COR2A-B view spans 2.5-15 R@. The temporal evolution of the top, middle and bottom panels is shown in 
movie2.mpg, movie3.mpg and movie4.mpg, respectively. 




Fig. 3. TB and PB images for STA (first row) and for STB (second row). They were taken at 18:22:20 UT. Pre-event 
TB and PB images were subtracted. Field of view spans from 2.5-15 Rq. 




Fig. 4. Schematic of the 3D models used in the modeling of the WL CME (red) and ray (blue). 
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Fig. 5. 3D geometrical modeling of the CME observations of STA/COR2 (upper panel) and ofSTB/COR2 (lower panel). 
The left column shows the TB images for each spacecraft, where a pre-event image was first subtracted. The middle 
column shows synthetic TB images of the CME from our 3D modeling whereas the right column shows overlays of the 
CME wireframe(red) on the corresponding TB images. Observations taken at 14:52:20 UT for STA and at 14:52:47 for 
STB. 
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Fig. 6. Schematic of the orientation (dashed magenta line) and of the angular width of the CME (solid magenta lines) 
as determined from the 3D modeling of Section 3. The plot also supplies the relative positions of STB, Earth and ST as 
well as the corresponding POS (blue, green and red solid lines respectively) during 9 April 2008. 
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Fig. 7. 3D geometrical modeling of the ray observations of STA/COR2 (upper panel) and of LASCO/C2 (lower 
panel). The left column shows TB images for each spacecraft, where a pre-event image was first subtracted. The middle 
column shows TB images of the ray from our 3D modeling whereas the right column shows overlays of the rectangular 
slab model wireframe(red) on the corresponding TB images. Times of the observation: 17:52:20 UT for STA/COR2 and 
17:50:04 for LASCO/C2. 




Fig. 8. Same as in Fig. [7] but for a ray model with twice the depth and the thickness of that of Fig. [7] 




Fig. 9. Determining excess masses(densities) along the ray (dashed line) from the COR2A observation at 17:52:20 UT. 
A pre-event frame taken at 09:22:00 UT was subtracted. 
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Fig. 10. Density along the ray as a function of heliocentric distance (squares). Density profile of the pre-CME corona on 
9 April 2008, 10:06 UT from LASCO C2 PB data inversion (dashes). Density distributions for a hydrostatic isothermal 
corona at 2 MK (dots). 




Fig. 11. Overlay of the CME (red cross) and ray (blue cross) central longitudes and latitudes determined from the 3D 
modeling of Section 3 on an MDI magnetogram taken during 1 April 2008, when the host AR (i.e., NOAA 1098) was close 
to the central meridian as viewed from SOHO. The dashed blue box has a length/height equal to the depth/ thickness of 
the WL ray as determined from our 3D modeling. Its orientation cannot not be determined from our modeling (see the 
discussion of Section 3). 
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Fig. 12. Combined TB view of the CME (green) and of the ray (orange). View from the ecliptic plane (left panel). View 
34 degrees above the ecliptic plane (right panel). Both views are from a distance of 100 Rq. 



